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The possibility of making an optically large (many wavelengths in diameter) object appear invisible has been 
a subject of many recent studies. Exact invisibility scenarios for large (relative to the wavelength) objects 
involve (met a) materials with superluminal phase velocity (refractive index less than unity) and/or magnetic 
response. We introduce a new approximation applicable to certain device geometries in the eikonal limit: 
piecewise-uniform scaling of the refractive index. This transformation preserves the ray trajectories, but leads 
to a uniform phase delay. We show how to take advantage of phase delays to achieve a limited (directional and 
wavelength-dependent) form of invisibility that does not require loss-ridden (meta)materials with superluminal 
phase velocities. © 2012 Optical Society of America 
OCIS codes: 000.0000, 999.9999. 



The attainability of optical invisibility devices (here- 
after, cloaks) has been a subject of active research and 
some debate in the recent years, immediately following 
the demonstration of a microwave cloak [1] based on 
the transformation optics [2] (TO) methodology. Recent 
progress in this area can be attributed to the develop- 
ment of electromagnetic (EM) metamaterials with un- 
usual properties such as strong magnetic response and 
large, controlled anisotropy across the spectrum, from 
radio [3,4] to terahertz [5] and even optical [6] frequen- 
cies. 

While TO is not the only methodology proposed for 
invisibility [7,8], so far it is the only known path to 
cloaking of very large (relative to free-space wavelength 
Ao) objects. To date, proposed and demonstrated TO 
cloaks based on passive, linear, reciprocal media require 
at least three [2, 9-12] EM properties that are hard to 
find in natural substances. One key ingredient is the large 
and precisely controlled anisotropy of the refractive in- 
dex [2,10] (RI). Implementations of TO recipes generally 
require anisotropic materials, with the exception of con- 
formal maps (CM). The usefulness of CMs for cloaking 
designs is fundamentally limited by several factors, in- 
cluding the preservation of the conformal modulus of the 
transformed domain, lack of rotational symmetry [11,12], 
and the absence of useful CMs in three dimensions. In 
2D, lacking radially symmetric CMs, the CM subset of 
TO can at best offer cloaking for one propagation di- 
rection, as was demonstrated by Leonhardt [13]. We re- 
fer to such devices as directional cloaks here. CM cloaks 
also suffer from other issues, such as the impossibility to 
conformally map a finite-area domain onto a smaller do- 
main, which leads to an inevitable impedance mismatch 
at the finite radius cutoff [11, 12]; such issues are easily 
avoided in non-conformal TO scenarios with anisotropic 



media. 

The second and more fundamental cloaking require- 
ment is superluminal phase velocity, or RI n < 1. Gen- 
erally, the requirement is n < no, where no is RI of im- 
mersion medium [14]; we consider only free space cloak- 
ing, no = 1. This requirement persists even in CM-based 
cloaking scenarios [11-13,15]. To implement n < 1, at 
least one of the principal values of the e and /i tensors of 
a non-bianisotropic medium would need to be less than 
unity. Transparent media with n < 1 are necessarily 
dispersive — otherwise they could be used to transmit 
information-carrying signals with superluminal speed, a 
clear contradiction with special relativistic causality. By 
virtue of Kramers-Kronig relations, one may deduce that 
effective media with e < 1 (or \i < 1) must be lossy be- 
cause of their dispersion. Photonic band gap media oper- 
ating beyond the effective medium regime in a high-order 
Bloch band have been proposed as superluminal velocity 
substrates for lossless, all-dielectric TO devices [16]; how- 
ever, they still suffer from the dispersion and impedance 
mismatch limitations. 

The third exotic property needed for TO invisi- 
bility designs is magnetic response. In 3D, a non- 
unity magnetic permeability (/i) is needed both for 
impedance matching as well as birefringence suppression 
in anisotropic media. In 2D, for in-plane wave propaga- 
tion, one may treat the decoupled TE and TM polar- 
izations separately, requiring the management of only 
a subset of the constitutive tensor elements. For TM 
waves, the magnetic permeability fi zz is a scalar. In the 
short- wavelength (eikonal) limit, one may further sim- 
plify the material properties of a TM-polarization cloak 
and thereby eliminate magnetic response [9,11]. Elimina- 
tion of non- unity \i from cloaking is valuable particularly 
for shorter wavelengths, where traditional magnetic ma- 
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terials have weak response and artificial magnetic meta- 
materials lead to a strongly dispersive and lossy fi [6]. 
Impedance mismatch can be eliminated or reduced us- 
ing an immersion medium with a properly chosen RI [14] ; 
however, this technique is not useful for practical appli- 
cations requiring cloaking in free space or low-density 
fluids. 

In this Letter, we concentrate on the simplified case of 
directional cloaks [13], further restricted to 2D in-plane 
propagation. Our main goal is to show that superluminal 
phase velocity can be eliminated using a variant of the 
eikonal approximation, namely, the scaling of RI. Addi- 
tionally, we show that the Fabry-Perot resonance (FPR) 
of the entire structure can be used to avoid the require- 
ment of non-unity \i from directional TM-wave cloaks. 




Fig. 1. (color online). The effect of RI scaling in a TM- 
polarization cloak with a = 1 m, b/a = 0.5 and B/a = 1: 
(a) Field plot (H z ) at A = (2a)/40; no RI scaling (fa = 
1). Black lines are EM flux streamlines, coincident with 
ray traces in the eikonal limit, (b) Same as (a) but with 
RI scaling factor fa = 2 (sufficient to make ei 5 2 > 1); 
A = (2a)/40. (c) Same as (b) (fa = 2), but at A = 
(2a)/34.2. (d) Measure of visibility, cfnf-, vs frequency in 
the index-scaled cloak with fa = 2; abscissa is Ac/)/ (2ir) = 
ko(2a)(fa — l)/(27r) (number of phase delay cycles per 
cloak length 2a). 

Our cloak geometry is shown in Fig. 1(a). The choice of 
geometry is not unique; the only requirement is that all 
exterior boundaries of the cloak are normal and parallel 
to the incidence direction. The 2D Helmholtz equation 
is modeled using the finite element method with COM- 
SOL RF solver. The diamond-shaped arbitrary- medium 
object of length 2a = 2 m and height 2b = 1 m cen- 
tered at the origin of Cartesian coordinates is coated by 
a perfectly conducting (PEC) sheet. The PEC-covered 
diamond is surrounded by a cloak with rectangular ex- 
terior boundary (length 2a and height 2B = 2a). 

The cloak consists of six homogeneous triangular 
pieces of material, labeled Region I, II, V — V" and IF 



in Fig. 1(a). Using the standard TO prescription [2], the 
material properties in regions I, II are given by 

f l =1 f l = f l = r f l = 1 + r 2 n l = 1 • 

xx " xy yx ' yy ' ' r*zz ' 

^xx X? ^xy ^yx ^, €yy 1/X> f^zz X> (X) 

where \ = B j(B — b) > 1 and r = b/a. Regions V — V" 
and IF are filled with a mirror-imaged material from 
the region I and II. These pieces implement a piecewise- 
linear, continuous transformation, 

I : x = .x, y' = y — rx — 6, 

II : X ' = X , y> = X (y - b ), (2) 

and likewise in regions V — V" and II', the result of which 
is the appearance of the diamond-shaped PEC object 
as a flat, zero-thickness, PEC sheet. In regions I and 
F, the transformation is area-preserving; thus, /j} zz = 1. 
The principal values of the dielectric tensors are e\ 2 = 

1 + ^ ± r^J 1 + 2 j-; note that e\ < 1. In regions II and 
IF, the transformation increases the area of the cloak 
inwards, thus allowing the object to shrink. Inevitably, 
n\\ > 1. In the same domains, the principal values of the 
dielectric tensors are e^ 2 = x ±l ] notably, eJ, 1 = l/x < 1- 

Next, we propose a scaling transformation that ensures 

i,ii ^ , 
e{ 2 > 1. 

To mitigate the need for RI or dielectric constant of 
less than unity, we introduce an approximation that pre- 
serves ray trajectories without necessarily maintaining 
the correct phase. This control is achieved by a uniform 
scaling of RI throughout the entire device, while simul- 
taneously preserving the wave impedance. The exterior 
index remains no = 1 in our procedure, which is differ- 
ent from the immersion approach [14]. This procedure 
creates two issues: (a) RI discontinuity on the surface 
of the device that causes unwanted reflections, and (b) 
phase discontinuity between the rays traveling inside and 
outside the cloak. 

To solve the issue (a), we have chosen our device geom- 
etry and function such that the RI discontinuity is per- 
pendicular to the wave incidence. At normal incidence, 
only the wave impedance affects the reflection coefficient, 
and RI jumps are invisible. This technique, of course, 
would not be applicable to omnidirectional cloaks, which 
are beyond the scope of this Letter. For the cloak shown 
in Fig. 1(a), the RI needs to be scaled up by at least a 
factor of fa = x = 2 to ensure e^ 1 > 1. 

Issue (b) causes interference between the rays passing 
close to the upper and lower boundaries of the cloak (at 
y = ±B). The interference along those boundaries can be 
suppressed by placing a pair of geometrically thin PEC 
sheets at y = ±5 spanning the entire cloak length, from 
x = —a to +a. The result is seen in Fig. 1(b). From one 
view angle, these PEC sheets are themselves invisible to 
TM waves in the eikonal limit, thereby preserving the 
cloak function. 

Even with the PEC shields, there is still significant in- 
terference past the cloak (Fig. 1(c)). This interference is 
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strongly suppressed at discrete frequencies, as seen from 
Fig. l(c,d). In Fig. 1(d), we plot a measure of visibil- 
ity denned as the norm of the magnetic field perturba- 
tion integrated over a line in the shadow of the cloak: 



(7NF 
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ds, where H z q is the amplitude of 



the incident wave. In our calculations, the line of inte- 
gration is placed at x = 1.1a, and its length is limited 
to 3B. The quantity cjnf has the units of length (in 
2D), and therefore can be regarded as a near-field scat- 
tering width. Evident from Fig. l(b,c,d), the minima of 
<jnf are good indicators of scattering suppression. One 
would expect that these minima should correspond to 
frequencies at which the phase delay between the inside 
and outside rays is an integer multiple of 27r; however 
the presence of diffracted waves causes the minima to be 
detuned from those positions. 





Fig. 2. (color online). Directional TM-polarization cloak 
without magnetic response (ji zz = 1): (a) Field plot (H z ) 
at v = 5.396 GHz. (b) Visibility figure of merit (ctnf) 
as a function of frequency, showing a series of FPR-like 
transmission resonances. 

The subluminal cloak with e > 1 (Fig. 1) still requires 
ji > 1. While a medium with \i > 1 can theoretically 
have arbitrarily small dispersion and negligible loss, such 
properties have only been found in materials at low fre- 
quencies. To tackle the magnetic response issue, we recall 
that in the geometric optics (GO) limit, the scaling of the 
wave impedance does not change ray trajectories. Arbi- 
trary scaling of the impedance profile was used by various 
authors and is known as the eikonal limit [1, 10, 11, 16] of 
TO. However, previous works considered only continuous 
impedance scaling, with the idea that, even in the GO 
limit, impedance discontinuities create reflections. Omis- 
sion of reflected rays on impedance-mismatched inter- 
faces led some authors to believe that CM-based cloaks 
may perform well at arbitrary incidence angles [12, 13]. 

While impedance discontinuities would increase visi- 
bility, we may still take advantage of the quasi-one di- 
mensional geometry of our structure. In ID, perfect in- 
visibility can be easily achieved — though only for a dis- 
crete set of Ao — using the well-known phenomenon of 
FPR. A lossless homogeneous slab exhibits 100% trans- 
mittance at certain Ao, regardless how large or small the 
impedance mismatch. Although a rectangular cloak is 
not exactly a ID object, we may expect that FPR-like 
resonances could strongly reduce the forward-scattering 



amplitude (i.e., the shadow). By virtue of the 2D optical 
theorem [17], then, the total extinction cross-section is 
also strongly suppressed. 

Full-wave simulations of the original (index-unsealed) 
cloak whose impedance was adjusted to have \i zz = 1 in 
all regions (Fig. 2) demonstrate that, indeed, numerous 
such resonances exist, and they strongly reduce visibility 
of the structure at select Ao- These resonances can be uti- 
lized as long as the length of the cloak 2a exceeds Ao/2. 
Their abundance in the GO limit, Xo/a —> 0, makes it 
easy to choose a close-to-desirable Xo/a ratio of the cloak. 
We emphasize that the use of either approximation pre- 
sented here creates a fundamental bandwidth limitation 
for the cloak — a trade-off to eliminating dispersive and 
lossy media. 

This work was supported through a Multidisciplinary 
University Research Initiative, sponsored by the U.S. 
Army Research Office (Grant No. W911NF-09-1-0539), 
and partially supported by the U.S. Navy through a sub- 
contract with Swamp Works (Contract No. N00167-11-P- 
0292). The authors are grateful to Nathan Landy (Duke 
University) for useful discussions on feasibility of direc- 
tional cloaks. 



References 

1. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, J. B. 
Pendry, A. F. Starr, and D. R. Smith, Science 314, 977 
(2006). 

2. J. Pendry, D. Schurig, and D. R. Smith, Science 312, 
1780 (2006). 

3. Y. Urzhumov, W. Chen, C. Bingham, W. Padilla, and 
D. R. Smith, Phys. Rev. B 85, 054430 (2012). 

4. D. Huang, Y. Urzhumov, D. R. Smith, K. H. Teo, and 
J. Zhang, Phys. Rev. B 111, 064902 (2012). 

5. H.-T. Chen, W. J. Padilla, J. M. O. Zide, A. C. Gossard, 
A. J. Taylor, and R. D. Averitt, Nature 444, 597 (2006). 

6. A. K. Sarychev and V. M. Shalaev, Electrodynamics of 
Metamaterials (Singapore: World Scientific Publishing, 
2007). 

7. A. Alu and N. Engheta, Phys. Rev. Lett. 100, 113901 
(2008). 

8. S. Tretyakov, P. Alitalo, O. Luukkonen, and C. Simovski, 
Phys. Rev. Lett. 103, 103905 (2009). 

9. W. Cai, U. K. Chettiar, A. V. Kildishev, V. M. Shalaev, 
and G. W. Milton, Appl. Phys. Lett. 91, 111105 (2007). 

10. Y. A. Urzhumov, F. Ghezzo, J. Hunt, and D. R. Smith, 
New J. Phys. 12, 073014 (2010). 

11. Y. A. Urzhumov, N. B. Kundtz, D. R. Smith, and J. B. 
Pendry, Journal of Optics 13, 024002 (2011). 

12. Y. Urzhumov, N. Landy, and D. R. Smith, J. Appl. Phys. 
Ill, 053105 (2012). 

13. U. Leonhardt, Science 312, 1777 (2006). 

14. B. Zhang, Y. Luo, X. Liu, and G. Barbastathis, Phys. 
Rev. Lett. 106, 033901 (2011). 

15. N. I. Landy, N. Kundtz, and D. R. Smith, Phys. Rev. 
Lett. 105, 193902 (2010). 

16. Y. A. Urzhumov and D. R. Smith, Phys. Rev. Lett. 105, 
163901 (2010). 

17. L. J. Boya and R. Murray, Phys. Rev. A 50, 4397 (1994). 



3 



Informational Fourth Page 
References 

1. D. Schurig, J. J. Mock, B. J. Justice, S. A. Cummer, 
J. B. Pendry, A. F. Starr, and D. R. Smith. Metamaterial 
electromagnetic cloak at microwave frequencies. Science, 
314:977, 2006. 

2. J. Pendry, D. Schurig, and D. R. Smith. Controlling 
electromagnetic fields. Science, 312:1780, 2006. 

3. Yaroslav Urzhumov, Wenchen Chen, Chris Bingham, 
Willie Padilla, and David R. Smith. Magnetic levita- 
tion of metamaterial bodies enhanced with magneto- 
static surface resonances. Phys. Rev. B, 85:054430, 2012. 

4. Da Huang, Yaroslav Urzhumov, David R. Smith, 
Koon Hoo Teo, and Jinyun Zhang. Magnetic superlens- 
enhanced inductive coupling for wireless power transfer. 
Phys. Rev. B, 111:064902, 2012. 

5. Hou-Tong Chen, Willie J. Padilla, Joshua M. O. Zide, 
Arthur C. Gossard, Antoinette J. Taylor, and Richard D. 
Averitt. High-speed terahertz modulation from metama- 
terials with embedded high electron mobility transistors. 
Nature, 444:597, 2006. 

6. A. K. Sarychev and V. M. Shalaev. Electrodynamics of 
Metamaterials. Singapore: World Scientific Publishing, 
2007. 

7. A. Alu and N. Engheta. Multifrequency optical invis- 
ibility cloak with layered plasmonic shells. Phys. Rev. 
Lett, 100:113901, 2008. 

8. S. Tretyakov, P. Alitalo, O. Luukkonen, and C. Simovski. 
Broadband electromagnetic cloaking of long cylindrical 
objects. Phys. Rev. Lett, 103:103905, 2009. 

9. W. Cai, U. K. Chettiar, A. V. Kildishev, V. M. Shalaev, 
and G. W. Milton. Nonmagnetic cloak with minimized 
scattering. Appl. Phys. Lett, 91:111105, 2007. 

10. Y. A. Urzhumov, F. Ghezzo, J. Hunt, and D. R. Smith. 
Acoustic cloaking transformations from attainable ma- 
terial properties. New J. Phys., 12:073014, 2010. 

11. Yaroslav A. Urzhumov, Nathan B. Kundtz, David R. 
Smith, and John B. Pendry. Cross-section comparisons 
of cloaks designed by transformation optical and opti- 
cal conformal mapping approaches. Journal of Optics, 
13(2):024002, 2011. 

12. Yaroslav Urzhumov, Nathan Landy, and David R. 
Smith. Isotropic-medium three-dimensional cloaks for 
acoustic and electromagnetic waves. J. Appl. Phys., 
111:053105, 2012. 

13. U. Leonhardt. Optical conformal mapping. Science, 
312(5781):1777-1780, 2006. 

14. Baile Zhang, Yuan Luo, Xiaogang Liu, and George Bar- 
bast at his. Macroscopic invisibility cloak for visible light. 
Phys. Rev. Lett, 106:033901, 2011. 

15. N. I. Landy, N. Kundtz, and D. R. Smith. Design- 
ing three-dimensional transformation optical media us- 
ing quasiconformal coordinate transformations. Phys. 
Rev. Lett, 105:193902, 2010. 

16. Yaroslav A. Urzhumov and David R. Smith. Transfor- 
mation optics with photonic band gap media. Phys. Rev. 
Lett, 105(16):163901, Oct 2010. 

17. Luis J. Boya and Robert Murray. Optical theorem in N 
dimensions. Phys. Rev. A, 50(5):4397, 1994. 



4 



